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I. INTRODUCTION 



General requirements. — From a physicist's standpoint, the most interesting part of quan- 
tum computing research may well be the possibility to probe the boundary between the 
quantum and the classical worlds. The more macroscopic are the structures involved, the 



better. So far, the most "macroscopic" qubit prototypes t 



"pes that hav 



lave been studied in the lab- 



oratory are certain kinds of superconducting qubits 12, 13|, 14, 15| . To get a feeling for how 
macroscopic these systems can be, the states of flux qubits which are brought in a quantum 
superposition corresponds to currents composed of as much as 10^ — 10^ electrons flowing in 
opposite directions in a superconducting loop. Non-superconducting qubits realized so far 
are all essentially mzcroscopic. 

The advantage of quantum superconducting circuits (QSC) as implementation of qubits 
is, first and foremost, due to the macroscopically quantum coherent ground state of super- 
conductors, which (a) supports non-dissipative flow of electrical current and (b) suppresses 
or outright eliminates low-energy elementary excitations. The latter property counterweighs 
the dangerous effects of a huge number degrees of freedom, which exist in the solid state and 
which would otherwise severely limit or totally destroy the quantum coherence necessary for 
the operation of a quantum computer. In this way the primary advantage of QSC allows 
the realization of the secondary one: as a solid state based device, a superconducting qubit 
and a more complex QSC can be more easily scaled up, can have significant density, and 
can be produced using a well developed set of design and fabrication methods. Moreover, 
since these circuits are macroscopic, it can be simpler to manipulate and read their state. 

The necessity to operate at low temperatures is not a disadvantage of QSC, since almost 
all the quantum information processing requires low temperatures in order to suppress the 
effects of noise. 

DiVincenzo criteria. — For the following, we will formulate an abbreviated set of require- 

il 

ments to a QSC serving as a register of qubits, condensed from the "DiVincenzo criteria" [6[ : 

1. A qubit is a two- level quantum system, which can be controlled and read out. 

2. A quantum computer is a set of qubits, where, in addition, certain two-qubit oper- 
ations can be performed. 

3. The decoherence time (i.e. the time during which the system maintains quantum co- 
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herence) must exceed the time necessary to perform a single- or two-qubit operation 
by a wide margin. 

A qubit state can thus be described by a vector in a 2-dimensional Hilbert space spanned by 
the eigenstates {|0), |1)} of a Pauh operator . Usually these states are chosen to coincide 
with distinct classical states of the system (e.g. corresponding to different electrical charge 
on a capacitor; different directions of current in a circuit) and which are easily distinguished 
in a measurement. "Controlled" means that an arbitrary unitary transformation U can be 
applied on an arbritary state of the system: \ip) — > U\il)). The "readout" requirement 
means that projection of on some directions in {|0), |1)} can be realised and that the 
resulting classical states can be distinguished. The state of a A^-qubit quantum computer 
belongs to a 2^-dimensional Hilbert space, (S)^i{|0)a) |l)a}) ^i^d the "two-qubit operations" 
demand that, for at least some qubit pairs a, 6, one could perform arbitrary unitary transfor- 
mations in the corresponding Hilbert subspaces, {|0)a, |l)a}®{|0)6, \ ^)b} in order to entangle 
these qubits. The last requirement, that of a long enough decoherence time, is necessary in 
order to preserve the quantum information through-out a computation. 

The major success of quantum information theory has been to show that quantum error 
correction [3] is possible and that therefore the computation time is not limited by the 
coherence time. As long as the error probability per gate is below a certain threshold, one 
is able to do arbitrary long quantum computations. In the literature, the most often cited 
threshold is 10~^, but this number is extremely sensitive to the assumptions about the noise 
and about the structure of the computer. It can be orders of magnitude smaller 8| or 
bigger , depending on the various assumptions that are made. Unfortunately, the number 
of 10~^ is widely cited, irrespective of these important assumptions. 

The DiVincenzo requirements 1. and 2. are clearly satisfied by a system whose "pseu- 
dospin" Hamiltonian is 

H = ^ha-aa + ^Sa-^ab-^b + (l) 

a a,b 

where = (c"a5 ^a) is the vector of Pauli operators acting on qubit a. In this expression, 
ha (the effective field acting on qubit a) and iab (the strength of qubit-qubit coupling) should 
be externally controllable to sufficient extent. 

Alternative approaches to quantum computing. — It is not our goal to discuss here the 
origin and the necessity of DiVincenzo criteria, for which one is better directed to handbooks 



on quantum computing (e.g. 
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Nevertheless it is useful to indicate certain approaches 



to quantum computing which allow to simplify these requirements, and would be especially 
well suited for an implementation using QSCs. 

The standard approach to quantum computing is based on quantum circuits, meaning 
the consecutive application of quantum gates (one- and two-qubit unitary operations) to 
the system, leading in the end to a highly entangled and fragile quantum state, where the 
solution to the problem is encoded. (A good illustration of the complexity of realizing a 



quantum circuit is given by the famous paper 



12| , reporting the factorization of the number 



15 using the room-temperature NMR approach. The operation required the application of 
around 300 vr-pulses and 30 7r/2-pulses.) The DiVincenzo criteria reflect the contradictory 
requirements to a quantum computer implementing the quantum circuits approach: on the 
one hand, the qubits must be maximally protected from the external influences to preserve 
quantum coherence; on the other hand, the need to do precise time-domain manipulations 
on them leads to complex control-and-readout circuitry, which will introduce decoherence 
into the system. 

The alternative schemes of quantum co.uputmg, =uch as cluster state Q, topological Q 
Il5l . Ilq . and ground state [17] quantum computing, avoid some of these difficulties. The 
most straightforward is adiahatic quantum computing (AQC) 18|, ll9(]. In AQC, the solution 
is encoded in the ground state of the system evolving under an adiabatically slow change of 
a control parameter. There are constructive proofs that for any standard quantum circuit 
one can realize the corresponding AQC with a local Hamiltonian and with no more than 
polynomial overhead 20|, l2l|] . Remaining in the ground state to a significant extent protects 
the system against relaxation and dephasing [22] . The precise time-dependent control over 
specific qubits is no longer an issue for AQC. Therefore a solid state-based, and in particular. 



QSC-based, implementation of an AQC is feasible 23 



24 



an exact realization of the prescriprions of Refs. 20|, l2l|] may prove difficult in practice, while 



25| . One should keep in mind that 



a deviation from it can lead either to exponential, instead of polynomial, overhead, or to 
the system ending up in one of the excited states. The question whether useful results can 



be still obtained in this case was discussed in 



26]. 



Josephson effect. — We have seen that a qubit must have two well defined levels that are 
used as logical states |0) and |1). However, in practice, very few systems in nature are defined 
by only two levels. To get around this, we can take two states in a nonlinear system and 
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treat them as an effective localized spin 1/2 (qubit). Due to nonlinearity, the transitions to 
other states can be made practically negligible. In electrical circuits, the natural solution is 
a superconducting Josephson junction 27], |28|: the only known nonlinear anc? nondissipative 
electrical element. 

Josephson junctions are formed by two superconductors separated by a weak link (a 
tunneling barrier, constriction, bridge, etc) (see Fig. [1]). "Weak" means that the probability 
amplitude for an electron to pass through the link (e.g. the tunneling matrix element K in 
case of a tunneling barrier) must be small enough to be considered as a perturbation. The 
coherent nondissipative current is carried by the Bose-condensate of Cooper pairs of electrons 
with (in most superconductors) opposite spins. A bulk superconductor is characterized by 
a position-dependent complex superconducting order parameter A which up to a factor can 
be written as y^exp(z</)). In this expression, rtg is the average density of electrons in the 
condensate, and is related to the velocity of supercurrent Vs oc V0. The supercurrent 
density is given by js = UsCVs- The order parameter plays the role of macroscopic wave 
function of the superconductor. 

Following the approach due to Feynman (see [s^]), valid in case of a low transparency bar- 
rier, we neglect the variation of the order parameter within each superconductor and describe 
the system by a two-component "wave function", \1/ = (y^ri7^exp(20i), y/n7^exp(202))^- 
Now, assume that there is a voltage diffence V between the superconductors, which means 
that the energy difference between them is 2eV (since the Cooper pair has charge 2e). Then 
the Schrodinger equation for \1/ can be written as 



(2) 



After a simple calculation we find that dug i/dt = {2Kns i/h) sin(02 ~ 4>i)', di'Pi ~ 4'i)/dt = 
2eV/h, that is, the celebrated formulae for the DC and AC Josephson effect: 

I J = h sin 0; = 2eV/h. (3) 

Here the phase difference = 02 — 0i. The first equation in ([3]) describes the non-dissipative, 
equilibrium, coherent flow of electric current through the barrier, determined only by the 
phase difference between the superconductors and the properties of the junction (all of the 
later being swept into the critical current Ic). In some cases, other functional dependences 
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than sin occur and these have interesting qubit apphcations that will be briefly discussed 
later on. The second equation ([3]) tells that if the phase difference is not constant, the current 
flow will be accompanied by a flnite voltage drop across the junction, and vice versa. 

The equilibrium current can be obtained from the appropriate thermodynamical potential 
of the system: Ij = dE/d^, where $ = ^Q(f)/2TT has the dimensionality of the magnetic 
flux, and $o = ^/2e is the superconducting flux quantum. Therefore we add to the energy 
of the system its Josephson energy, 

Ej{(f)) = [ d^Ij{ct>) = -^cos0 = -Ejcos0. (4) 

J ^TT 

From the AC Josephson effect and according to the relation V = LI, we see that the 
Josephson junction can be considered as a nonlinear inductance Lj: 

" (5) 



2e/c cos 

The fact that this inductance can be tuned (and even change sign) by the stationary phase 
difference across the junction is very useful for the applications. 

As an example of this quantum inductance, consider the rf-SQUID (Fig. [2b), i.e., a 
superconducting loop of inductance L interrupted by a single Josephson junction. The 
quantization condition , which follows from the single- valuedness of the order parameter, 
connects the phase drop across the junction to the total magnetic flux through the loop via 

+ 27r(<l>,.-L/j)/$o = 27rn. (6) 

Here is the total phase drop across all Josephson junctions in the loop (if there are more 
than one), is the external flux through the loop, and L is the self-inductance of the loop. 
The response of the current in the loop to small changes in and therefore the effective 
inductance of the SQUID, depends on the value of ^x- 

The total energy of a Josephson junction also includes its electrostatic energy: 

E = QV2C + Ej{(j)) = EciNc - N^f + E^). (7) 

Here Q is the electric charge on the junction's capacitance, C , and Ec = 4e^/2C the charging 
energy (the Coulomb energy of one Cooper pair). The charge Q is expressed through Nc and 
Ng. The former is the number of Cooper pairs having tunneled through the junction while 
the latter depends on the system's conflguration and its electrostatic interaction with its 



surroundings. E.g., when there is a voltage biased electrode of voltage Vg and of capacitance 
Cg to the junction (see Fig. [2^), we have A^^^ = CgVg/2e + 6Ng where ^A^^^ is the effective 
charge on the junction due to any other sources (like charged impurities) in the proximity. 
While Nq is a discrete number, Ng can take arbitrary values. 

To certain extent, iVj^ and can be considered as conjugate variables, similar to momen- 



tum and position 28|. We can therefore quantize Eq.([7]) by replacing the classical variables 



(j) and Nc with operators satisfying the commutation relation 

[Nc,^] = -t; Nc = -tdf, (8) 

we omit the hats on the operators further on. The situation is in reality more complex than 
that, because Nc cannot take negative values, and is only defined modulo 2tt (see e.g. 



Ref. 29|, §1.4.3). Nevertheless for large enough systems with Nc ^ 1, Eq. ([8]) holds. Given 



this, the junction Hamiltonian becomes 

H = Ec{-td^-Ngf-Ejcos(l). (9) 

Taking Ng = and expanding the cosine around = 0, we obtain the Hamiltonian of a 
linear oscillator with the eigenfrequency c^o: —^^o4' = (p = {ih)^'^[[(f), H], H], that is, hwo = 



\j2EcEj. In practice, this plasma frequency is in the range of tens of GHz. 

II. BASIC TYPES OF SUPERCONDUCTING QUBITS 

Three main classes of superconducting qubits, all based on Josephson junctions, have 
been theoretically studied and tested experimentally. These are phase, flux and charge 
qubits. The relation between the parameters hujQ, Ec and Ej and the way these qubits are 
biased distinguishes between the different types. 



Phase quhit. — Phase qubits [30|, |3l|] operate in the phase regime, which is defined by 
Ej ^ Ec- In this situation, the Josephson term dominates the Hamiltonian Q. These 
qubits consist of a single current-biased Josephson junction (Fig. EJi). The Hamiltonian of 
the system can be written as 

H = -Ecdl - Ej cos - ^0 = -Ecdl - Ej (^COS + y^0^ , (lO) 

where the last term describes the effect of the external bias current This is the Hamil- 
tonian of a quantum particle in a tilted washboard potential, see Fig. [2]i. Since Ej ^ Ec, 



in the absence of bias, h = 0, each energy well contains many almost equidistant quantized 
levels, and the tunneling between different minima (0 = 0, ±2tt, ±47r, . . . ) is negligible. On 
the other hand, the Josephson junction cannot support a nondissipative current exceeding 
If.. Indeed, when 1^ > I^., the potential in flTU]) no longer has local minima. The resulting 
solution with nonzero (0) is called resistive state: the current is accompanied by finite volt- 
age [see Eq. ([3])], and is dissipative due to the normal resistance of the junction (which we 
neglected so far). 

In the subcritical regime (in practice, when If, is about 0.95/c — 0.98/c ) there remain only 
a few quantized levels in every local minimum. The tunneling out of the first two levels in 
a given potential through the barrier is still small enough, and these can be taken as qubit 



states |0) and |1) [30|, l3l|- When ~ Ic the phase ~ 7r/2; expanding the Josephson 



potential in ffTOj) to the third order in = — 7r/2, we find 

Uj{^) = Ej [{1 - h/h)^ - <^V6] . (11) 

The plasma frequency of the resulting cubic oscillator (i.e. the frequency of small os- 
cillations in the quadratic potential approximating f|TT]) near its minimum) is ujp{Ih) = 
ujq{1 — {Ih/ IcYY^^. The transition frequency between the two qubit states is cuqi ~ 0.95ujp 
and depends on the bias current. The effective qubit Hamiltonian ( fTOl) becomes sim- 
ply H = —^^^a^, corresponding to a spin-1/2 in a field along the z direction. The x- 
and/or y- components, necessary for unitary rotations of the qubit, can be implemented by 
adding an oscillating component to the bias current at the qubit transition frequency Uqi. 



h = h,dc + -^bac "^os ^01^ + Kac^'^'^^oii- The effective Hamiltonian becomes [32| 



f^Ol z , / ^ Kmc X , / ^ K-ac 



2 V 2ujmC 2 V 2ujQiC 2 ^ ' 

The DC component of the bias current sets tuoi, while the cosine and sine quadratures 
produce the rotations around the x and the y axis, respectively. The nonlinearity of the 
potential (|TT|) is crucial, otherwise the AC current pulses would also excite transitions be- 
tween higher levels, and out of the Hilbert space spanned by |0) and In other words, 
the system would no longer be a qubit. 

Advantages of phase qubits include the simplicity of design, good control over the param- 
eters in the Hamiltonian, simple readout, weak sensitivity to charge and flux noise in the 
system, and scalabihty. One of the relative disadvantages is the sensitivity to low-frequency 
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(i.e. 1//) noise of critical and bias current. This noise can be produced, for example, by 
microscopic defects in the insulators of tunneling barriers. Another is the absence of the 
optimal point, i.e., the degeneracy point where symmetry would protect the qubit from 
certain kinds of noise and, as we will see for other types of qubits, significantly increase its 
coherence time. This is what has so far limited the system to short coherence times (T2 ~ 80 



ns 



33i). 

In phase qubits cr^'^-operations must be realized with AC pulses. This has the disadvan- 
tage of requiring complex high-frequency low-noise circuitry down to the very low operating 
temperatures (~ 10 mK). On the other hand, using AC rather than DC pulses allows for 
protection by low-pass filtration of the qubit from the low frequency (1//) noise coming 
from the bias lines that can be significant. 

Flux qubit. — Another possibility to realize a qubit in the phase limit Ej ^ Ec is by using 



a degeneracy between two current-carrying states of an RF-SQUID (see Fig. [Jb) j3J]. In 
writing the Hamiltonian of the system, the only difference with respect to the Hamiltonian 
of Eq. is that we now must take into account of the magnetic energy, but can neglect 
the effect of off-set charges Ng-. 

H = -Ecdl - Ejcos<f) + El{<P - 4>xf jl. (13) 

Here (^x = 27r$a;/^o is the reduced flux through the loop, and the inductive energy scale is 
given by El = ^I/AttL. The situation differs from the previous case because now, if (px ~ tt, 
the potential energy formed by the last two terms of Eq. f|T3l) has two almost degenerate 
minima. As illustrated in Figure [2)d, these states correspond to persistent current in the 
loop circulating in opposite directions, and are conveniently used as the |0) and |1) states 
of a qubit. 

Tunneling between the two potential wells is enabled by the charge term in Eq. ( fT3i) . As 
a result, in the subspace {|0), |1)} the effective qubit Hamiltonian is 

H = -'-a^ + ^a^, (14) 

where e oc (l|i/|l) — (0|//|0) is the energy bias between the two states, which is tuned by the 
external flux; A oc {0\H\1) is the tunneling amplitude. When e ^ A, the eigenstates of ([H]) 
coincide with |0) and |1). At the degeneracy point, e = ( ^^.^ = tt ), these eigenstates are 
"bonding/antibonding" states |±) = (|0) ± \l))/\/2. The expectation value of the current 



circulating in the SQUID loop, /, is the same in |+) and |— ): (+|/|+) = (— ). As a 
result, these two states cannot be distinguished by observing the current, or the magnetic 
field it produces. Therefore no external degrees of freedom that couple to magnetic field can 
interfere with the qubit when it is operated at this optimal point; the coherence time of the 
qubit can thus be significantly enhanced. Therefore it is convenient to always work at this 
special operating point and use the corresponding eigenstates as effective logical states of 
the qubit: \0l) = |+), = |— ). While the Hamiltonian (HM enables all qubit rotations 
with DC pulses only, transitions between these states can also be realized by applying RF 
flux pulses at the transition frequency, same as in a phase qubit. 

Note that no difference in observed current means that it is impossible to readout the 
state of the qubit by looking at current or magnetic fleld produced by the loop. One can 
either move away from the optimal point for the readout, as was done in earlier experiments, 
or to use dispersive readout discussed later. 

The major problem with the RF-SQUID design is that in order to have a potential energy 
landscape with two well-defined minima, both the inductance of the SQUID loop and the 
Josephson energy must be large. First, this suppresses the tunneling; in the experiment 34 1 
it was impossible to produce coherent superpositions between the lowest states in the wells, 
and instead the higher states, close to the top of the barrier, had to be used; transitions 
to these states were induced by resonant RF pulses. Second, due to large inductance the 
states |0) and |1) produce quite large magnetic fields (these two states differ by about 
one flux quantum $o) and the minute deviations from the optimal point lead to strong 
coupling to the environment and fast dephasing. This problem is solved in the flux qubit 
design of Ref. [35|, l36|. Fig. [Jb, where the RF SQUID is replaced by a loop of negligible 
inductance L but with three Josephson junctions. In this case the total flux through the loop 
practically coincides with the external flux, and the flux quantization condition ([6]) becomes 
01 + 02 + 03 + 27r$2;/$o = 27rn (here the phase drop across junction i is (pi). It leaves only 
two independent phases (e.g., 0i and 02)- The Josephson energy of the system provides a 
2D potential landscape, -E'j(0i, 02; 0z) = —Eji cos0i — Ej2 cos 02 — Ejs cos(0s + 0i + 02), 
plotted in Fig. [Jb and containing suitable double well structures. The two qubit states 
again produce circulating currents in the qubit loop. However, the generated flux is now 
much smaller, and the two states are separated by only ~ 10^^$o- As a result, unwanted 
coupling to the environment is much weaker, with appropriately longer coherence times. The 
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tunneling barrier can be also lowered, enabling tunneling between the lowest states in each 
minimum. The Hamiltonian of this qubit can be also cast in the form f|T^ . and the same 
optimal point considerations apply. Note that a consistent derivation of the limit L ^ 
from Eq. flT^ is rather involved 37 1. 

While the fluxes produced by the persistent currents in |0) and |1) only differ by about 
10~^$0) the difference of currents themselves still involve ~ 10^ — 10^ single-electron states. A 
state of the form (|0) ± |1))/ thus involves a superposition of millions of charges flowing in 
two directions at once. Superconducting flux qubits are truly macroscopic quantum systems. 
(A more detailed analysis of how to define what is a macroscopic quantum system is given 



m 



38 



39|.^ 



Advantages of flux qubits are their weak sensitivity to charge noise, comparatively easy 
readout, the ability to work at or near the degeneracy point (optimal point). Moreover, 
in principle, all operations on the qubit can be performed using DC pulses only. Their 
main disadvantage is a very strong dependence on the junction parameters. The tunneling 



splitting A is proportional to \^EjEc exp[— ^y Ej / Ec], while Ej itself depends exponentially 
on the barrier thickness. Such fluctuations as 1// noise in the critical current will thus have 
a big effect on the qubit. 

Charge qubit. — The charge qubit (Fig. [2^) operates in the opposite regime from the 
previous two types of qubits, the charge limit Ec ^ Ej. Now the states with definite charge, 
differing by a single Cooper pair, are the working states of the qubit, and the Josephson 
term is the perturbation. Historically, this the first superconducting qubit to have been 
experimentally implemented 40, |4l|]. This is not completely surprising given that, in a way, 
it is the "most microscopic" of them all. 



In order for a single Cooper pair to make a difference, two conditions must be met |28 |. 
The charging energy Ec ~ e^/C must far exceed both the thermal energy, ksT, and the 
linewidth of the charge states, AE ~ h/RC, due to their finite lifetime determined by 
the effective resistance and capacitance of the system. Otherwise, the charge states are 
smeared out and cannot form a good basis for a qubit. The second condition leads to 



28l . §7.1, this becomes a problem. 



the requirement i? > ^/e ~ 6 kil. As noted in Ref. 
because at GHz frequencies dictated by small capacitances C, the tunneling resistance will 
be effectively shunted by small impedance ~ 50 — 100 Q of the connecting leads (that is, bulk 
superconductors forming the banks of the Josephson junction). The charge states will be 



11 



washed out by the quantum particle number fluctuations. As a result, it is easier to observe 
charge effects in a small superconducting island separated by high-impedance tunnel barrier 
from the bulk superconductor, as shown in Fig. [2^. The Hamiltonian for a charge qubit can 
be written as the sum of the electrostatic energy and of the Josephson energy: 

H = Ec{Nc - Ngf - Ejcos0. (15) 

In the absence of Josephson energy, the states \Nc = 0) and \Nc = 1) are degenerate when 
Ng = CgVg/2e = l/2. Higher-energy states can be neglected; then on the Hilbert subspace 
spanned by these two states, the operator Nc takes the form Nc = — cr^). The Josephson 
term lifts the charge degeneracy; in this basis it can be written as Eja^ /2. The Hamiltonian 
again takes the pseudospin form 

H = ~Ecil-2Ng)a^-^a^, (16) 

allowing all unitary rotations with DC pulses only. At the charge degeneracy point, A^^ = 
1/2, the eigenstates of this Hamiltonian are coherent superpositions of states (|0) ± |l))/-\/2 
differing by a single Cooper pair. Like in the flux case, it is best to operate at this optimal 
point, and it is convenient to make a transformation to the basis of logical qubit states 
(|0) ± \l))/V2: (T^ ^ -a' and cx^ ^ cr^ . Then, denoting Ng = 1/2 + Ng^acCos{ujt), we have 

H = + EcNg^ac cos(u;t)cT^. (17) 

Now the Josephson energy plays the role of transition frequency for the qubit and AC voltage 
Ng^ac can be used to induce transitions between these states. 

The advantages of the charge qubit operating at the optimal point are similar to those of 
the flux qubit. The former however does not have the exponential dependance on parameters 
that the flux qubit has. One of its main disadvantages is very strong sensitivity to charge 
noise. This can be mitigated to some extent by working in the intermediate regime, Ej < Ec, 



which was realized in the so-called "quantronium" 42| and, more recently, the "transmon" 



qubit 



43|. 



III. COUPLING QUBITS 

Coupling by passive elements. — By passive elements, we refer to the circuit elements with 
resonance frequency much higher than the |0) to |1) transition frequency of the qubits they 
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are coupling. As a result, the coupling element will remain in its ground state at all times. In 
this category, the most commonly studied (and experimentally exploited) coupling elements 
are the usual linear elements of electrical circuits: capacitances and inductances. Simplest 
passive couplings for phase (a), charge (b) and flux (c) qubits are shown in Fig. [31 For phase 
qubits, joining one electrode of each junction by a capacitance yields an interaction term 



44| 



H... = ^J^^^:^^,y,l (18) 

between the qubits a and h. Here ujQi^j is the transition frequency of qubit j. If the difference 
between the qubit eigenfrequencies, cuoi.a — ^m,b, is much bigger than the coupling strength 
h^ujQi^a^Qi^hC X / C , then this interaction term acts as a perturbation on the qubits and can, 
for all practical purposes, be ignored. Indeed, by moving to the interaction representation 
with respect to the Hamiltonian of noninteracting qubits, we will see that this term is 
fast oscillating and averages to zero, and is therefore negligible (so called Rotating Wave 
Approximation, RWA). On the other hand, if t^oi^a = ^m,b) the term (ITS]) cannot be neglected 
and leads to the qubit-qubit coupling. Therefore the coupling can be turned on and off by 
taking the qubits in and out of resonance with each other (for phase qubits, by tuning the 
qubit bias currents). If we use only DC controls, there is an important difference between 
the diagonal couplings (cr^a^) and off-diagonal ones {cf^ctI and (t'^^ctI). While the latter can 
be turned on and off by tuning the qubits, the former is always on. 

In practice, one cannot neglect finite inductances of connecting elements (roughly L ~ 
1 nH per 1 yum). Therefore even in Fig. [3^ the qubits are coupled through a harmonic 



oscillator of frequency c^osc = \/^I^Cx- Taking the capacitance to be 1 fF 33] and L =1 
nH, its frequency is ~ 100 GHz, about 10-20 times larger than the typical qubit transition 
frequencies. Therefore we have indeed a passive coupling and can safely neglect the finite 
inductance. In one experiment on capacitively coupled phase qubits with large ~ 5 pF 



the effects of the excited LC mode were observed 



45|. 



The capacitive coupling of two charge qubits. Fig. [3]d, leads to a similar coupling term 



e2 



^int = l^fa<. (19) 

to which the same considerations apply. 

Finally, Figure [3t shows two flux qubits coupled through a mutual inductance M. Since 
external flux leads to a change in energy bias ea^ /2 in the single qubit Hamiltonian Eq. (fT4l) . 
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this coupling is 



(20) 



where Ij is the circulating current in qubit j. In the eigenbasis of the qubits at the optimal 
point, this also takes the form 



One can use a bona fide LC circuit as a passive coupler [46|, |47|, |48| (like in FigJlt)); 
now several qubits can be connected to a single circuit {quantum bus), and only those in 
resonance with each other will couple. If the qubits themselves cannot be tuned, one needs 



more complex passive coupling elements (e.g. 



49 



;unable sign was recently realized experimentally 51 



50[). For flux qubits, coupling with 
521]. The coupling scheme used in 



5l| is shown in Fig. Ht. The external flux through the coupler loop, $c = fc^o, changes 



the phase shifts across the junctions shared between the coupler and the qubits a and b. 
The resulting shift in the Josephson energy of the whole system leads to the (diagonal) 
qubit-qubit coupling Jcr^a^ with 

n I'Uc) 



C) 



2e /2 - /(/c 



(21) 



Here Ic is the critical current of the shared junctions, and I{fc), I'ifc) are the circulating 
current in the coupler loop and its derivative with respect to the coupler flux resp.; 
are the persistent currents in the qubits. The coupling sign changes with the isgn of the 
derivative I'{fc), and the experimental dots in Fig. are in a good agreement with the 
formula. 

Use of AC controls opens still more possibilities. The diagonal coupling can be effectively 
cancelled by the techniques of refocusing borrowed from NMR j53| . The off-diagonal coupling 
can be tuned with RF pulses while staying at the optimal point. The general idea is to use the 
frequency control. For example, detune the qubits from each other, so that 6 = tuoi.a — ^oi,b 
exceeds the coupling strength; the coupling is effectively off. RF irradiation of each qubits 
at the difference frequency 6 will effectively turn on the interaction 



55 



56 



57|, 



58 



. Several different 



59 



60| 



frequency-control schemes have been recently proposed 

Coupling by active elements. — By these we understand the circuit elements whose typical 
transition frequency are comparable to those of the qubits they couple, and which therefore 
undergo transitions to excited states. In other words, now the coupling between the qubits 
is realized through their exchange of real, rather than virtual, excitations with the coupler. 
Examples of couplers which can be used in the active regime are shown in FigureHl In Fig.H^, 
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charge qubits are coupled through a current-biased Josephson junction (CBJJ) [4J,|6l|. Like 
in the case of a phase qubit, the junction can be considered a tunable nonlinear oscillator, 
with only the two lowest states involved, and is coupled to either qubit through gjO'ab'^osc- 
The coupling strength Qj depends on the coupling capacitance Cxj between the qubit and 
the oscillator. Since this interaction is off-diagonal, when the oscillator is strongly detuned 
from the qubits, the oscilator-qubit interaction is turned off. By sequentially tuning the 
oscillator in and out of resonance with the two qubits, it is possible to entangle these qubits 
while leaving the oscillator unentangled [4J]. In this way, several qubits can be coupled to 
the same oscillator, which plays the role of the data bus. Since it is the coupling bus and 
not the qubits that are tuned, the qubits can remain at their optimal point throughout the 
manipulations. The same approach works for other types of qubits as well (Fig. |Ud). 
For a more exhaustive review of coupling schemes the reader is directed to Ref. 62 1. 
Qubits coupled to oscillators: circuit QED. — The natural way superconducting qubits 
and oscillators couple to each other led to the insight that superconducting quantum circuits 
allow to realize an analogy to quantum optics in general 63|], and to strong coupling cavity 
quantum electrodynamics (cavity QED) [6J, l65| in particular. This was pointed ou t by 



several authors 



44 



66 



67 



68 



69|. In this section, we will focus on the work of Ref. 70|. 



As shown in Fig. [5], this proposal is based on a charge qubit (green) fabricated inside a 
high quality superconducting transmission line resonator (blue). In the lumped element 
description, the transmission line resonator can be seen as a series of LC circuits connected 
to external leads by capacitances Cq. The finite length / of the resonator sets up its char- 
acteristic frequency. The electromagnetic field transmission between the input and output 
ports of the resonator is maximal at this frequency and its harmonics; with / = 25 mm, 
these resonances are in the microwave range. 

Close to one of these resonances, the resonator Hamiltonian is = huj^a'^a, where 
ujr = l/"\/ZC*, and a (a^) are the annihilation (creation) operators. The qubit is biased not 
only by a DC or AC current (that can be applied on the input port of the resonator), but also 
biased by the voltage across the LC circuit. This voltage can be written as Vic = Kms('^^+'^)) 
where V^^^ = ^/hjIlJ2C is the rms value of the volta ge i n the ground state. In practice, this 
rms voltage can be as large as a few /iV 43|, |7l|, |72|, |73|] . Inserting this in the Hamiltonian 



of a charge qubit DC-biased at the optimal point f[T71) . we find H 



-a'- 



where g = eV^^JJ g / Cy, is the qubit-oscillator coupling strength, with Cg the oscillator-qubit 
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coupling capacitance and Ce total capacitance of the qubit. Writing cr^ = o""*" + cr and 
neglecting fast oscillating terms (RWA), we can write the qubit +oscillator Hamiltonian as 

H = hujra)a + h^(j^ -ng{a)(j- + (j+a), (22) 



with fkjja = Ej. This expression is the Jaynes-Cummings Hamiltonian of quantum optics 63 1 



which is traditionally used to describe the physis of a two- 



evel atom of frequency Ua coupled 



to the photon field inside a cavity of frequency Ur 6J, l65|]. Due to this correspondence one 

calls this architecture circuit QED. 

Circuit QED has several advantages over cavity QED. First, due to the small volume of 

the transmission line resonator and to the large effective dipole moment of the charge qubit, 

the coupling strength g can be made much bigger. Second, the positions of real atoms inside 

a 3D cavity fluctuate in space, leading to a fluctuations in the coupling g, which interfere 

with quantum information processing. In circuit QED, the position of the qubit in the 

resonator is fixed. Third, it is possible to fabricate several qubits inside a single resonator 

I I 

and to entangle these qubits by using the resonator as a quantum bus i74i]. Circuit QED 
is thus interesting both for quantum information processing applications and as a tool to 
study new regimes of quantum optics. 



IV. READOUT OF SUPERCONDUCTING QUBITS 

To be useful as qubits, the quantum states of QSCs must allow an efficient readout. In 
practice, the readout circuitry will be microfabricated in a proximity to the qubit. The 
challenge is thus to have this circuitry strongly coupled to the qubit during measurements 
(to have a fast measurement on the scale of the qubit 's Ti) but strongly decoupled during 
computations (to prevent additional decoherence and relaxation of the qubit). 

For phase qubits, the readout circuitry is an integral part of the qubit design. The 
measurement is based on the large difference between the tunneling rate Fj of the logical 
states \i) 3^. As illustrated in Fig. Eli, to realize a measurement, the DC bias current is 
rapidly increased in such a way that while levels |0) and |1) remain in the well, the level 
|1) becomes very close to the top of the potential barrier, so that Fi ^ Fq. As a result, 
if the qubit is in the state it will rapidly decay in the continuum [75], producing the 
resistive state and thus a voltage of the order of the superconducting gap (~ 1 mV) across 
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the junction. Detection of this voltage measures the qubit state. An important drawback 
of this is the generation of quasiparticles due to the finite vohage. These will affect the 
qubit itself and lower the coherence time of neighboring qubits. Therefore, in more recent 
experiments the junction forming the qubit was embedded in a loop of finite inductance. 
Now the system tunnels from the local potential minimum not into the continuum of states, 
but in a large but finite potential well. As a result, there is no voltage drop across the 
junction, but a change of flux in the qubit loop by $o- The latter is detected by a DC- 
SQUID [76j]. Remarkably, if the qubit is not kept biased long enough for the state |1) to 
decay with certainty, a partial measurement is realized 77|] , accompanied by the probability 
amplitude concentrating in the state |0). 

For charge and flux qubits, the first readout methods to be experimentally realized de- 
tected directly the charge (flux) associated with the qubit states. For charge qubits |41| . 
this was first realized by capacitively coupling a superconducting single-electron transistor 
(SET) [40|], a highly sensitive charge meter, to the island of the qubit. The radio-frequency 
version of the superconducting SET (RF-SET) TSj has also been used in more recent ex- 



periments 



79 



80 



81| . By working at RF rather than at DC, the RF-SET can be much 



faster than regular SETs and is also less sensitive to 1// noise due to charge fluctuations. 
For flux qubits, the direct readout was realized by inductively coupling a DC-SQUID to the 



qubit 



34 



36|. 



A problem with the above mentioned approaches is that they do not work at the optimal 
point, where, as discussed earlier, the observable values of currents (fluxes) or charges in 
both qubit states coincide. Away from the optimal point the system is exposed to extra 
decoherence and relaxation, especially if the qubit 's environment is highly structured 



76 



. A 



partial solution is provided by the phase-charge readout developed by the Saclay group |42|. 
A different approach, dispersive readout, is based on the observation that, while the currents 
(charges) at the optimal point are the same in both qubit states, their derivatives with respect 
to the control parameters are not. This can be done by coupling the qubit to a resonator and 
measuring the latter 's susceptibility, dependent on the quantum state of a qubit. In case of a 
relatively low-Q lumped LC circuit far away from resonance (strong detuning) this approach 



was developed and implemented for flux qubits by the Jena- Vancouver co^ 



82 



83 



84| (with direct inductive coupling) and by the Delft group 



85 



86 



laboration 



51 



87| (with the 



coupling through a DC SQUID). In a similar way, measurement of the quantum capacitance 
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was also realized with lumped resonators in the very strong detuning regime with charge 



qubits 



88, 



89|. 



In the remainder of this section, we will discuss in more detail t 
circuit QED, where this approach is especially effective 



431, 



71 



le dispersive readout in 



73| . It works by going 

to the dispersive regime, where the qubit-resonator detuning A = h{uOa — t^r) = is much 
larger than the coupling strength g. In this situation, the Jaynes-Cummings Hamiltonian of 
Eq. fl22l) can be replaced, at the optimal point {Ng = 1/2) and to second order in the small 
parameter g/A, by the effective Hamiltonian 



(23) 



where Ua = oJa + X ^"^^ ^X = 9^ I The term x^^^^^'^ dZSj) acts as a shift of the cavity 
frequency that depends on the state of the qubit. Indeed, in this situation the effective 
resonator frequency is no longer lOr, but is now ± X depending on (cr^) = ±1. Using this 
frequency pull, it is possible to read out the state by measuring the phase and amplitude 
of an RF signal transmitted between the ports of the resonator. These quantities indeed 
contain information about the state of the qubit {toI . 90 |. 

Dispersive readout has several important advantages over the approaches that directly 
measure charge or flux. First, the qubit can be operated anc? measured at the optimal point. 
This is because it measures not the charge on the qubit's island but rather the resonator 
frequency pull which is related to the quantum capacitance of the qubit [t^. Second, it is 
actually at the optimal point that this quantum capacitance (or quantum inductance in the 
flux case) is maximally different for the two states of the qubit. Thus, with this dispersive 
measurement, there is maximal signal at the optimal point. Third, there is no energy 
exchange between the qubit and the resonator systems during readout, and the resonator is 
measured in its eigenstate. This leads to what is known as quantum non- demolition (QND) 



measurement 



63| and is ideally suited for quantum information processing. Last, but not 



least, the readout is only active when signal is sent at the input port of the resonator, which 
realises the desired coupling/decoupling requirement for a quantum readout. 
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V. FURTHER READING 



The short format of this review allows us neither to include all the important topics con- 
cerning superconducting quantum circuits, nor to cover the ones included in due detail. We 
try to remedy this by providing here a list of additional sources, which should be accessible 
based on what have been discussed here so far. 



Reviews on SQC. — A concise and accessible description of different types o: 
ducting qubits is provided in Refill and, with a varied degree of generality, in 



471 . Recent pedagogical reviews are Refs. 



supercon- 



2 

62, 911. The 



The first comprehensive review is Ref. 
paper 9l| describes all qubits as examples of a nonlinear quantum circuit, with an emphasis 



on the phase qubit. The review 



62| contains also a description of abstract qubit manipula- 



tions (Sect. IV) and of the theoretical formalism necessary to analyze classical an d q uantum 



92 



9^ 



superconducting circuits (Sect. V and VI). Recent advances are reviewed in 

Superconductivity. — Introductory information on theory of superconductivity and 
Josephson effect can be found in [27|. The classic book [28| devotes a whole chapter (Ch.7) 



to the effects in sma^ 
is also covered in 



1 Josephson junctions highly relevant for the qubit theory. This topic 



29|, Ch.4. 



Theoretical formalism. -Detailed explanation of the formalism allowing a uniform descrip- 



tion o 
in 



94 



arbitrary quantum superconducting circuits, based on a Lagrangian approach, is given 
951]. The formalism is linked to the graph description of complex circuits, standard 



in electrical engineering. 



coupling are derived in [37|, l96| . 



vffective Hamiltonians for qubit manipulations and qubit-qubit 



Decoherence and dissipation. — Sour c es of decoherence in var i ous superconductin g qubits 



are discussed in ha. 



97 



98 



99 



lOG 



10 ll . 



102 



103 



104 



105 



106 



107 



108 



109 |. 



Recent experimental progress. — After significant improvements to the initial design, two- 



qubit quantum operations (equivalent to a univ ersal ViSWAP gate) were demonstrated in 



two capaciti yely 



coupled phase qubits [33|, lHOj. In particular, in Ref. 



33l . quantum state 



tomography lll| was realized to extract the density matrix representing the state of the 
coupled qubits. An experiment with two phase qubits capacitively coupled through a series 
LC resonator showed the evidence for entanglement between the qubits and the resonator 



mode 



451]. In single phase qubits the behaviour of quantum two-level systems (TLS), which 



are beheved to be the main source of 1// noise, was investigated, and coherent beats between 
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a TLS and a qubit were observed [Ty, llOq . This led to a proposal [48| to use such TLSs as 



naturally formed qubits, with the phase qubit playing the role of a quantum bus. 

With flux qubits, much effort was devoted to the development of structures suitable 



for adiabatic quantum computing. After the initial proposals 23|, |2J], coupling between 



two 



112], three 25| and four 8J] flux qubits was realized in the quantum regime, and 



entangled eigenstat es w ere formed. We already mentioned tunable coupling between two 



flux qubits 



51 



52 



llSj. Time- domain tunable coupli ng b etween flux qubits and a simple 



quantum protocol were experimentally realized in R ef. 114| . Multiphoton Rabi oscillations 
were observed in a flux qubit 1151 ]. In Ref. 116| a flux qubit was measured using the 
Andreev probe, that is, measuring the electrical conductance in a non-superconducting (in 
this particular case, silver) wire connected by superconducting leads to two points of the 
qubit loop. The normal current through the wire is carried by so-called Andreev levels in a 
normal conductor sandwiched between superconductors (see |29|, §A.2). The conductance 
is therefore sensitive to the phase difference, and allows to determine the state of the qubit 
with a minimal backaction. 

The longest relaxation and dephasing times for charge qubits (8 /is and 500 ns, respec- 
tively) were measured with a qubit fabricated inside a transmission line resonator (circuit 
QED) [t^. Recently, the same design was used to non-destructively measure the ph oton 
number population of a resonator [43]. In flux qubits dephasing times as long as 4 /xs 117 ] 
were observed. 

A new latching readout based on the non-linear dynamics of a large Josephson junction 
was implemented 118| . Distinguishing features of this approach are that it is dispersive 
and that, once the readout as occurred, information about the logical state of the qubit is 
encoded in the dynamical state of the junction for a long time. This allows, in principle, 
to accumulate enough signal about the state of the qubit to overwhelm any noise in the 



circuitry. This type of latching readout was also applied to flux qubits [STt^ 
Exotic superconducting qubits. — "Phase slip" qubit proposed in Ref. 



119| replaces the 



tunneling Josephson junction of an rf SQUID qubit by a thin superconducting wire (bridge), 
which plays the role of the tunnel barrier. The geometric inductance of the loop can be kept 
small, because now the main role is played by the kinetic inductance which relates the 
supercurrent increase to th e change in the superconducting phase difference, L^Ig 



"Andreev level" qubit 



271 



1201 ] replaces the tunneling Josephson junction by another type of 
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weak link, a quantum point contact, i.e., a constriction, in which only a few conducting 
modes can fit (see Ch.4). The Andreev levels are formed inside the constriction and 

carry the Josephson current. In either case, the effective tunneling matrix element should 
not exponentially depend on the microscopic properties of the wire (point contact). It is 
therefore expected that such qubits would be easier to fabricate, and they would be less 
sensitive to 1// n oise, than a flu x qubit. "d-wave qubits" using high-Tc superconductors 



were proposed in 



121 



122 



123| . Due to their (i-wave pairing symmetry, the Josephson 



junctions with such superconductors may have an intrinsically degenerate ground state and 
therefore do not require external circuitry to stay at the optimal point. The gr ound state 
degeneracy in YBCO Josephson junctions was experimentally verified in Ref. 12J|. The 
feared decoherence from nodal quasiparticles could have been overestimated Recent 
observations o : 



and BiSCCO 



macroscopic quantum tunneling in ci-wave junctions based on Y BCO 



126 



127(1, as well as coherent transitions in an YBCO junction 



125| 



128) . lend 



credibility to these conclusions and make the possibility of high-Tc qubits more realistic. 
Interesti ngly , YBCO junctions with vr-phase shifts were also recently incorporated in RSFQ 
devices 1291], which, of course, work in a classical regime. Due to the interplay of magnetic 



and superconducting ordering, a Josep 
an equilibrium phase difference tt 96 



ison junction with a ferromagnetic barrier can have 



ISOj. Inserting such a junction in a flux qubit loop 



(a "SFS qubit") is equivalent to putting in it a $o/2-flux quantum and thus produces the 
intrinsic degeneracy of the ground. 

In conclusion, we can confidently state that superconducting qubits and more general 
SQCs have a serious potential as "quantum hardware"; but it is their unique suitability for 
probing the boundaries between quantum and classical behaviour of nature, which makes 
them so attractive objects of research. 
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FIG. 1: Josephson junction and Josephson effect. The state of the jth superconductor is charac- 
terized by the "superconducting electron density", n^j-, and superconducting phase, In case 
of tunnehng Josephson junction, shown here, both the critical current Ic of the junction and its 
effective capacitance C depend on thickness of the insulating layer, through which the Cooper pairs 
tunnel. 
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FIG. 2: (a) Simplest charge qubit [411]. (b) RF SQUID (single-junction flux qubit [34']) and its 
potential energy, when the external magnetic flux through the loop = ^o/2. The flux quantum 
$0 = h/2e. (c) Three-junction flux qubit 3^ and its Josephson potential profile as a function 
of two independent phase differences (see text) when = Dots mark degenerate states, 
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and its Josephson energy ("washboard 



between which tunneling occurs, (d) Phase qubit 
potential"). Two metastable states indicated can be used as qubit states ]0) and ]1); Fq and 
Fi ^ Fq are the tunneling rates out of these states into the continuum. 
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FIG. 3: (a) Phase qubits and (b) charge qubits coupled through a capacitor, (c) Flux qubits 
coupled through the mutual inductance M. 
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FIG. 4: (a) Charge qubits coupled through a tunable bus circuit (red) (adapted from Ref. 



(b) Flux qubits coupled through an LC circuit (red), (c) Tunable coupling of two flux qubits 



(adapted from Ref. 



5 lip. The coupling is tuned between ferro- and antiferromagnetic by changing 



the magnetic flux, <I> = fc^Oi through the coupler (red). Dots show the experimental data. 
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FIG. 5: Circuit QED (adapted from Ref.Q): 

a charge qubit coupled to a strip line, and the 

simplified scheme of the system (inset). 
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